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Abstract 

In a recent paper S. Friedland and the author presented a formal ex- 
pression for Xd(p) of the monomer-dimer problem in d dimensions involv- 
ing a power series in p. We there presented the result of computations for 
the terms in the power series through the sixth power. I herein present 
the result for the seventh power term, in d — 2. An interesting feature 
of the new computation is that it would have been impossible without a 
new algorithm for computation of the Tutte polynomial, by Bjorklund, 
Husfeldt, Kaski, and Koivisto. 

In [?] the author presented a formal expansion for of the dimer problem 
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working with Shmuel Friedland this work was extended to the monomer-dimer 
problem, yielding a formal expansion 

MP) ~ ~(ph(2d) - P ln(p) - 2(1 - p) ln(l -p) -p) + (2) 

fe=i 

with 

ci(p) = \p 2 (3) 

c 2 (p) = (2p 3 + 3p 4 )/96 (4) 

c 3 (p) = (-5/ + 12p 5 + 8p 6 )/192 (5) 

We have come to believe that it is better to organize the expansion in @ as a 
power series in p. We thus write 

^ oo 

A d (p) ~ 2(pln(2d) -PHP) - 2(1 - P) In(l -p) -p) + ^ a k (d)p k (6) 

k=2 
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with 



1 1 



(7) 



03 (p) 



8d 
1 1 



(8) 



48 d 2 



<h(p) 



05 (p) 



1_ 1_ 5 1 
11 39 1 



(10) 



(9) 



a 6 (d) 



16 d 3 640 d 4 
11 11 



19 1 
1920 



(11) 



24 d 3 32 d 4 



Each of the three expansions is computed as far as it can be from knowledge 
of the kernels Ji computed to construct series ([1]). These Ji are given in (31)- 
(36) of [?], and repeated in (5.25)-(5.30) of [?]. The Ji were obtained by lengthy 
computations, by machine for i = 5,6. Once the Ji are known each of the three 
series above is given by an easy computation (a short run of a short Maple 
program). 

When I first constructed these series I viewed them as asymptotic, series ((TJ) 
valid in the limit d — > oo. And series (|6]) for p small enough. But we now believe 
series © converges for all physical values! I.e. d = 1, 2, • • • , < p < 1. As was 
pointed out in [?], (6)— (11) is true for d — 1; and for d = 2, (6)— (11) is strongly 
consistent with a numerical study of A2 (p) ■ 

For d = 2, the expression from (6)— (11), and the computation of the present 
paper, is ( first Ji = yfffg and therefrom ) 



The computation of the p 7 term involved 1000+ hours' running of a Maple 
program on my p.c. The computation of the contribution of each perturbation 
term was done in integer arithmetic, but to avoid ridiculous fractions at interme- 
diate steps, the sum of contributions was done in floating point. The interesting 
feature of the computation involved the computation of the expression ijj' c , see 
eq. (24) of [?] or eq. (5.21) of [?]. The overlap pattern of the dimers determines 
a graph with 7 vertices each representing a dimer, and edges between verticies 
whose dimers overlap. ip' c as defined in [?], when computed from its definition, 
requires 2 21 computer subroutines, 2 21 being the number of possible subsets of 
the edges. I learned from Richard Kenyon that i/j' c i s ; U P to sign, Tq(1,0), Tq 
the Tutte polynomial of the graph. And in [?] an algorithm is given which has 
an exponential in the number of vertices, instead of the number of edges, de- 



A2 (p) ~ ^ (P ln ( 4 ) - P Inp - 2(1 - p) ln(l - p) - p) 
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termining its computation time (to compute T"g)0 Use of this algorithm made 
the computation possible. 
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1 In fact I implement the "baseline algorithm" of [?], Section 3.1 therein, but using the 
basic zcta and Moebius transforms rather than their fast versions 
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